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Abstract 

Integrated forms of the one-loop evolution equations are given for the Yukawa 
couplings in the MSSM, valid for any value of tan (3, generalizable to virtually any 
number of Yukawa fermions, and including all gauge couplings. These forms turn out 
to have nice mathematical convergence properties which we prove, and we determine 
the ensuing convergence criteria. Furthermore, they allow to write down general 
sufficient and necessary conditions to avoid singularities in the evolution of the 
Yukawa couplings over physically relevant energy ranges. We also comment briefly 
on the possible use of these features for physics issues and give a short numerical 
illustration. 



1 Introduction 



Large Yukawa couplings often play an important role in spontaneous electroweak sym- 
metry breaking models such as the standard model, its extensions, and also in some 
alternatives to it. The existence of Infrared (IR) attractive fixed points |L| and effective 
IR fixed points in such a regime can be a benchmark in selecting the phenomenologi- 
cally viable theories in view of the measured mass of the top quark. Such considerations 
clearly favoured the minimal supersymmetric extension of the standard model (MSSM) 
H|, with its upper bound on the top mass, m t < (190 — 2WGeV) sin f3 J| over alternatives 
such as ||. 

In supersymmetry there is of course much more to it. More specifically, in the minimal 
supersymmetric extension of the standard model (MSSM) |J, the spontaneous breaking 
of the electroweak symmetry is generically driven by the large top quark mass in associ- 
ation with a soft supersymmetry breaking sector ||. Such a possibility can be looked at 
as the phenomenological facet of a [still to be discovered] deep connection between the 
origin of supersymmetry breaking and that of the electroweak symmetry breaking. Mean- 
while, it helps in correlating theoretically the very many free parameters of the MSSM, 
leading to quantitative estimates of the supersymmetric partners spectrum, of prime im- 
portance in guiding the experimental search for supersymmetry, and relating (at least 
qualitatively) the physics over as many orders of magnitudes as between the GUT scales 
and the electroweak scale. 

Of course, the key point in all the above issues is the way the various parameters 
"run", as dictated by the renormalization group equations. The RGE's are, in general, 
complicated coupled differential equations already at the one-loop level and one usually 
resorts to numerical methods to solve them 0. However, analytical solutions would be 
desirable for several reasons [besides the obvious one of allowing a better control over the 
structure of the running.] The radiative breaking of SU C (3) x SUl(2) x C/y(l) and the 
structure of the vacuum are controlled by the effective potential (EP) of the theory. The 
theoretical improvement of the functional form of the EP needs the analytical form of the 
running of involved quantities such as masses, couplings and fields which enter the game. 
If the analytical form of the EP beyond the tree-level were known exactly, or at least in 
an RG improved form beyond the naive loop corrections, one could hope to determine 
the relations the initial values of the parameters should fulfill, to break the electro-weak 
symmetry at the right energy scale and avoid in the same time color and charge breaking 
vacuaQ Another important application of the analytical solutions is the determination of 
the effective behaviour in the low energy regime (infra-red regions) independently (but 
within a domain) of the actual values in the ultra-violet, such as the the top Yukawa 
coupling effective fixed point [g], or the triviality bounds on the Higgs mass. 

Analytical solutions of the complete set of RGE's in the MSSM were known to one- 
loop order for small tan (3, strictly speaking in the case all Yukawa couplings are put to 
zero except for the top quark. Actually the structure of the coupled equations is such 

Mn the absence of such a knowledge, due partly (in the RG improvement program) to the existence of 
many different mass scales, one relies on rough approximations hoping that they encompass the leading 
behaviour ||. 



2 



that a necessary condition to solve them entirely is to be able to solve first for the Yukawa 
couplings. This is of course not enough, and initially one assumed also universality of the 
soft susy breaking terms at the GUT scale in order to solve for those quantities too. This 
assumption can however be relaxed [|Hj but still for small tan/?. It is thus natural to 
try to find exact solutions for the Yukawa sector for any value of tan /3, for comparable 
top and bottom Yukawa couplings, and also bringing in the game the r lepton Yukawa 
coupling as well in order to cope with the case of b — r unification [ IT| . Some attempts 
have been made (for instance in |12[) to solve generally the top-bottom Yukawa system 
which lead to implicit solutions provided one neglects the U{1)y gauge coupling (see also 



In the present paper we study some properties of the runnings of the Yukawa couplings 
as dictated by the RGE's to one-loop. The first aim is to provide suitable expressions for 
the exact solutions, which we call "integrated forms" . Although these expressions do not 
appear in closed forms, they are especially convenient because, if one insists on making 
them explicit, they come out as continued integrated fractions the convergence of which 
can be kept under control. Although our main results are a priori valid for any gauge 
theory with an arbitrarily extended Yukawa sector ( including the special case of the stan- 
dard model), we restrict most of the discussions and further illustrations to the case of the 
MSSM. We will give integrated forms, valid for any number of Yukawa couplings, of the 
general explicit solutions corresponding to the coupled renormalization group equations 
which read in the MSSM |TJ 



Gauge Couplings [gi with i — 1, 2, 3 and n g the generation number]: 

10 



-17 = T^hgf with 6 a 



dt 32tt 2 ^ 1 _i 3 

Yukawa Couplings [i = 1,2,3 generations]: 



1 - -rrn g , b 2 = 5 - 2n n , h 



9-2n n 





Y % 

u 


dt 


32tt 2 


dYj 




dt 


32vr 2 


dY{ 


n 


dt 


32vr 2 



k=gen \ / 

k=gen \ y 6 / 



mf+ e {(^ fc ) 2 +3(^) 2 }-3(^+^ 

k=gen 



;i-2) 

1.3) 

;i-4) 



Here the evolution parameter t is defined by t = Log(My/Q 2 ) where Mjj denotes 
some initial scale. Note that since a gauge coupling unification condition is not es- 
sential in the present study, we write the RGE equations in terms of the low energy 
SU(3) C x SU(2)l x U(l)y gauge couplings, respectively #3, #2 and g\. Note also that we 
assume here, and throughout the paper, flavour conserving (diagonal) Yukawa matrices. 



The rest of the paper is organized as follows. In section 2 we recall the known solution 
for large top Yukawa coupling and give the integrated form of the general solution valid 
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for any value of the top and bottom Yukawas. We then generalize those integrated forms 
to any number of Yukawa couplings, in particular to the top-bottom-r system. In section 
3 we give a proof of the convergence of these forms in both top-bottom and top-bottom- 
r cases. Section 4 is devoted to the question of avoiding Landau poles in the Yukawa 
runnings. There we give a generalization to the top-bottom case of some well known 
bounds, and establish necessary and sufficient conditions. Preliminary applications and 
comments are made in section 5 and conclusions and an outlook are given in section 6. 
An appendix contains some detailed proofs and technical material. 



2 Integrated form of the Yukawa coupling RGE's 

2.1 Large top quark Yukawa solutions: a reminder 

We are interested here in eqs. ( |1.2| - |1.4| ). They can be treated independently of the rest 
of the system, especially from the gauge couplings for which the running is determined 
a priori via eq. [This is no more true at two-loop order where the gauge and 

Yukawa equations become highly interwound.] When all Yukawa couplings except Y t are 
neglected, eqs. (|i~3l ) and ( |1.4|) become trivial while eq. (|1.2|) becomes of the Bernoulli type 
in the variable yt = Y 4 2 

j t Vt = fi{t)y t + by 2 t (2.1) 



where 



/iW = T^(t^ + 3 ^ 2 + T^ ' b = -^-2 (2-2) 

l07T z 6 9 l07T z 

and is easily solved to give 

Vt[t) - 1 — by Jq E(t')dt' (2 - 3) 

where 

E(t) = e /o/i(*')*' and y o = y t 2 (t = 0) (2.4) 



2.2 top-bottom case 

In the more general case where both Y t and Y b are kept in the game, but neglecting all 
other Yukawa couplings, eqs. (|I^ , |l~3| ) read, after the change of variables y t = Yf, y b = Y 6 2 , 

-j t Vt = h{t)Vt + ay b y t + by\ , ^y b = f 2 (t)y b + ay b y t + by\ (2.5) 
where fi(t) and b are given in eqs. (|2.2|) and 



1 46 o . „ o + 7 
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Ut) = -^(-fl, 2 + 3o 5 2 + -o 2 ) , a = (2.6) 

16tt 2V 3 ya y2 9 yiJ 16tt 2 V ; 



As far as we know, the system eqs. ( |2.5| ) is not treated in standard text books, and 
although it looks at first sight simple, we could not find a systematic way of relating it to 
a standard form Q. It is also relatively easy to solve the system up to first order in Y b in 
the region Y t 3> Yj>. This is already an improvement of the known solutions with Yj, ~ 0. 
It extends the numerical validity much further than tan/5 ~ 10. 

More importantly, this approximate solution gives a valuable hint for the structure of a 
suitable integrated form which can then be found by sheer guess and will be written down 
below. But first, the approximate solution can be obtained in the form y t (t) = y t (t) +S(t) 
where y~t{t) is given by eq. Q2.3| ) and linearizing when necessary the equations, in the regime 
yb(t), \5(t) \ <C 1. One then finds for y b 



1 - by° b J* E 21 (t')dt> 



Vb(t) = , o , t „ ( 2 - 7 ) 



where 



EAt) ~ ^Mjmw 1 (2 - 8) 

Ei(t) = efiWW i = 1,2 (2.9) 



and a slightly more complicated expression for y t given in appendix B. A little thinking 
then leads to the following form of exact "solution" we were looking for: 



V?E 12 (t) 
l-by?JoE 12 (t')dt> 

l-byUlE 2l {t>)dt> 



where 



Ei(t) 

'l-by° b JoE 21 (t')dt'Y/ b 
E 2 (t) 

:i-by?foEi2(t>)dtr /b 



Eu(t) = -, „ t , ; — r^r (2-12) 



E2i(t) = r, rTTT^TI^E (2-13) 



and y^ = Y t 2 (t = 0),y b = Y b 2 (t = 0) are any initial conditions. The reader can easily 
check that the solutions ( |2 . 1 0| , pTTT|) satisfy exactly Eqs.( |2.5| ) without any restriction or 

2 The situation would be much simpler if fi(t) = /^(t), in which case the equations can be solved 
by quadrature after some change of variables, leading though only to implicit solutions involving some 
hypergeometric functions O] 
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assumption about the magnitudes of the Yukawa couplings, i.e. for any value of tan/3. 
They resemble formally eq. ( [2.3|) of which they are a generalization. Of course, although 
our solutions, yt,Vb are now explicit in terms of E\ 2 and E 2 i, the latter are given only 
implicitly by Eqs. Q2.1^ , [2.13| ), which appear as coupled nonlinear integral equations. 
Therefore, the procedure is useful only if it provides us with a systematic (and hopefully 
quick) way to solve these equations within a given accuracy. It will be shown in sect. 3 that 
mere iterations achieve this goal. In fact, such iterations correspond to the truncations 
of the "continued integrated fractions" which naturally emerge as formal solutions of 
Eqs.([g£T3]), e.g.: 



E {t) = m 

< J " bM ^irhcrt )a,b 

- H J ° )a,b 

v-wj? iiUdt )a,b 



(2.14) 



2.3 Arbitrary number of Yukawa fermions 



In fact the above solutions are easily generalized to include any number of leptons and 
quarks. For instance, if one includes in the game a third Yukawa coupling, then Eqs.( |1.2| , 
|1.3| , |1.4| ) take the following form: 

-nVi = fi(t)yi + a n yf + a 12 yiy 2 + a 13 yxy 3 

^2/2 = /2(*)*/2 + «222/2 + 0-212/22/1 + 0-232/22/3 
372/3 = /3 0O2/3 + ^332/3 + «3l2/32/l + 0322/32/2 



dt 



The exact solution reads: 



(2-15) 
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(2.16) 



where 111,112 and M3 are defined through the implicit system 
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(2.17) 



and / Uj stands for f dt'uj(t'). 

In the interesting case of top-bottom-r system with y t = y 1: y b = y 2 and y T = y 3 one 
has in the MSSM 

6 4 
an — a 22 — — - - 2 ; a 33 



167T 2 167T 2 
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/»(*) = 1^(0? + ; ^(t) = e J?^ 



f2.i8) 



and fi, 2 {t), Ei 2 {t) as previously. It is interesting to note that in this case w r and u t 
are directly related via 

^ (2.19) 



-E3 ^■£'1 

a reflection of the fact that, in the MSSM, the running of yt and y T at one-loop order are 
mutually affected only indirectly through the running of y b) ( 031 = a i3 = 0) at variance 
with the non supersymmetric Standard Model (SM) case. 

Finally, the extension to more than three Yukawa couplings will not play any role in 
the present paper. Nevertheless, we give it for the sake of completeness in appendix C. 

3 Proof of convergence 
3.1 The top-bottom case: 

In this section we make a mathematical digression to study some useful properties of our 
solutions. Even though equations ( |2.12| , |2.13| ) give Eij only implicitly, they enjoy the 



property of defining a contraction mapping. This is about all what one needs to give 
a rigorous proof for the existence and uniqueness of the E^, and thus of the existence 
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and uniqueness of the solutions given in Eq.( [2.10 , 2.11|) . This proof will also be of prac 



tical use. It provides us with a criterion for the convergence of the truncated forms of 
Eq.( |2.14) ) towards the exact solution, and the rate of this convergence can be controlled 
so that a very good approximation will be obtained with a few (or even just one) iterations. 

For the sake of completeness, we recall here in simple terms the conditions required 
for a contraction mapping, and then prove that they are indeed satisfied in our case. Let 
us define 

U x {t) = (3.1) 



E 1 (t) 



and think of U\ and [/ 2 as forming a vector 




U(t) = rrTA (3-3) 



in some space St where the evolution parameter t remains in the interval < t < T for 
a given value of T. Then Eqs. (|2.12| , |2.13|) restrict the U-s to the positive region Ui(t) > 



(provided one stays far from the Landau poles), and, furthermore, define a mapping in 
this region, A: U \—>U' through 



U[{t) 



by° 2 JoE2(t')U 2 (t')dtr /b 
1 



(3.4) 

The idea now is to show that the mapping A shrinks uniformly, at each iteration, the 
"distance" between any two vectors in £ T (subject to the condition Ui(t) > 0). More 
precisely we will prove that there exists a positive constant number Kt < 1 such that the 
following inequality is satisfied 



|| U'-V \\< K T || U-V || (3.5) 

for any pair of vectors (U, V) belonging to St and satisfying [/,, Vi > ( i — 1, 2). Here 
|| . || is defined by 

|| U ||= max{ sup | U\{t) |, sup | t/2^) |} (3.6) 

0<t<T 0<t<T 
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Then, according to the "contraction mapping principle" , the existence of the (unique) 
solution of Eqs. (p.l2| , |2.13| ) in St is guaranteed, and the n th iteration A n U approaches 



this solution at least as fast as Kj,. 



To prove Eq. Q3.5|) one writes the following sequence of inequalities 

yU t dt'E j {t')\U j {t')-V j {t')\ 



um-v>(t)\ < 



[16tt 2 + 6y° min{/ 4 dt' E 3 (t')U 3 {t') , J* dt'E 3 (t')V 3 (t')}]V 6 

t 

< -f^ sup I U 3 {r) - V 3 {r) I / dt'E 3 {t') 

107T Z n< r <T Jo 



0<t<T JO 
/" _ _^ rT 



< i!h\\U-V\\ I dt'E 3 (t>) 
\bn z Jo 

(3.7) 

valid for i ^ j with i,j = 1,2 with y 1 = y t , y 2 = yb and where we plugged the actual 
values of the coefficients a, b in Eq. fl2.5p . The first inequality in ( |3.7| ) is derived from 
Eqs.( p.4|) and from the inequality, 

l< C jZ ■ OTWH ( 3 - 8 ) 



lv l + a 7 v l + /3 7 (l + min{a,/?}) c+r 
valid for any a,/3, larger than —1 and c > (here c = 1/6) 

the second from the positivity of V^, £/j and and the third from the definition 
( |3.6| ). Eq. (|3.7p immediately leads to Eq.(|3.5|) with 

K T = -r^maxiy" [ dtE 1 {t),y° b [ dtE 2 (t)} (3.9) 
\bn z Jo Jo 

The convergence condition Kt < 1 is easily met even when T is large enough to 
encompass the whole evolution range from the GUT scale to the Mz scale. For instance, 
if T « 66 and a^ T ~ 25, one needs Y®,Y b ° < 0(tt) where the Y® are the GUT scale 
values of the Yukawa couplings, to ensure convergence. These conditions are naturally met 
within the perturbative regime. We should stress, however, that these are only sufficient 
conditions. 



3.2 The top-bottom-r case: 



As we said previously the exact solutions ( |2.16| , p,17| ) to the generalized equations with 
three Yukawa couplings or more, is a direct generalization of eqs. (|2.10 , |2.11| ). Similarly 
the proof of convergence goes essentially along the same lines as in the previous section, 
when supplemented with the inequality 
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1 1 

< 



;i + ai y (i + (5 x y (i + a 2 y (i + (3 2 y 

a I ai — oii I £> I /?i — /?2 I 



l + minja!,^})^ 1 (1 + min^, /? 2 }) fe+1 

(3.10) 



valid for any a,b,ai,j3i > 0. Defining a mapping through eqs. fl2.17f ), one finds the 
convergence criterion 

Kt = ■ max .mo^'C - ^? / ^' ~ a ' lky ° k / E ^ < 1 ( 3 " 1:L ) 

(ijk)ctrc.perm.oj (12o) JO JO 

In the top-bottom-r case it reads 

1 T T T T 

Kt = j^ m a x{ % ° jf £ 2 ; y t ° jf ^ + y° jT E 3 ! 3^° jT E 2 } = 

j 

— maxK / + y° / E 3 ; 3y° b [ E 2 } < 1 

07T^ JO JO JO 
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(3.12) 



We see that the sufficient convergence criterion can become more severe in this case by 
about a factor 3 in the regime y° ~ ~ 



4 Avoiding Landau poles 

One question that can be clearly answered with the knowledge of analytical solutions is 
how to determine the conditions which guarantee that the values of the Yukawa couplings 
at some low-energy scale, say the electroweak scale, remain consistent with a Landau 
"pole" free theory up to a given high energy scale, typically a grand unification scale. 
The answer would be trivial if one starts from the high energy scale Yukawa coupling and 
runs down. Indeed in this case, it is clear from the general form of the solutions, eqs. (|2.10 
, |2~TT1 ) and the fact that b < 0,y°(= Y° 2 ) > and E{t)ij > 0, that one does not hit a 
Landau pole all the way down below the initial scaleQ The situation is more complicated 
if one starts from some Yukawa coupling values at a low scale and tries to run upwards to 
determine the corresponding values at a GUT scale. This is a phenomenologically typical 
situation if a model-independent reconstruction of the fundamental parameters is to be 
carried out, starting from the experimental input. 



3 Note that here we are only interested in poles which can occur explicitly in the Yukawa couplings. 
In practice one stays anyway far from gauge coupling Landau poles, given the running range relevant to 
our discussion. 



10 



Moreover since in the vicinity of such Landau poles the Yukawa couplings become 
very large, the conditions for avoiding these poles correspond in some cases to effectively 
attractive fixed points, such as the celebrated relation 



top 



;i90- 200GW)sin/? (4.1) 



in the MSSM valid when all Yukawa couplings are neglected in comparison to the that 
of the top quark [[|. When the Yukawa couplings are of comparable values such a cor- 
respondence becomes more involved, as can be seen for instance from the dependence on 
y® and y® in (|2.14j) [See also the discussion in section 5.2]. 

To illustrate the case, we start first with the solution for small tan (3, i.e. Y t ° 3> Y b ° ~ 0, 
given in eq.( [2.3|) . Writing it in the form 

E(t) + by t (t)tiE(t>)dt> y -> 
one sees immediately that y t (t) should satisfy 

m < -ij^h (43) 

for any value of t > in order to be consistent with the positivity of y° and ytif). Ensuring 
this positivity avoids automatically the Landau pole ( we rely here on the fact that b < 0, 
see eq.(2.2) .) Thus contrary to the initial value at a high scale t — 0, values of y t can 
not be arbitrarily chosen at lower scales, the maximum allowed in this case being simply 
given by eq. ( |4.3|) . With this in mind, it is straightforward to eliminate completely the 
dependence on y° in the running to get 

, M y t (to)E(t;t ) 

1 - by t {t ) J t0 E(t';t )dt' 

where 

E(t;t ) = ^- (4.5) 
Here yt(t ) is any initial value satisfying eq.( |4.3| ) for t = t . 
In the more general case, when y b is not neglected, one can also write, 

yb{to)E 2 i{t;t ) 



Vb(t) 
Vt{t) 



l-by b (t )flE 21 (t>;to)dt> 

yt{to)E 12 {t;t Q ) 
l-by t (t )flE 12 (t>;t )dt> 

(4.6) 
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with 



Eu(t', to) 



Ei(t;t ) 



(l-by b (to)flE 21 (t>;to)dtr /b 
E 2 {t;t ) 



(l-by t (to)f t t E 12 (t>;to)dt>y/ b 



(4.7) 



where 



E 1 (t;t )=e-l t ° dt '^ 



Ej{t) 



(4.8) 



Ejito) 



Note that here to > t corresponds to an initial energy scale lower than the running scale 
corresponding to t. In this case, however, it is no more possible to determine easily a 
sufficient and necessary condition (if any) on yb(to) and yt(to) to avoid the Landau pole. 
Indeed the sufficient and necessary conditions for a non singular running in the interval 
[T, t ] read 



as easily seen from eqs. (|4.6| ). However, E 12 and E 21 depend themselves on y® and y® so 
that the above conditions are highly implicit in y® and y®. 

Instead, one can immediately determine some necessary conditions, and, with some 
extra work, also some sufficient ones. We list these conditions below and refer the reader 
to the appendix for the detailed proofs. 

The necessary conditions: one can write a tower of pair of inequalities, each pair being 
a necessary condition weaker than the subsequent one, 




(4.10) 



(4.9) 
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u s „ glfa) 



Vb(to) < 



EM 
\b\$dtE 2 (t) 



(4.11) 



Vt{to) < 



Ei(t ) 



dtxExfa 



151 r*o 

'■' r (l-\b\y b (to)j£E 2 (t 2 ;t )dt 2 )"/» 

E 2 (t ) 



Vb(to) < 



\b\I, 



to 



dt x E 2 {t x ) 



(4.12) 



yt(t ) < 



Vb(to) < 



1 



\b\$dtE 12 (t;t ) 
1 

\b\ffidtE 21 {t;t ) 



(4.13) 



The limit of this tower of inequalities (assuming of course that the iteration converges) 



is precisely the necessary and sufficient condition eqs. ([4.9| , |4.10|) , as can be seen from 
eqs.( |4.7|) when written in a form similar to eq. (|2.14|) . 



The sufficient conditions: 



yt(t ) < 



Eiito) 



c(l + l/cY+ c \b\ ffidtE^t) 



(4.14) 



1 E 2 {to) 

mito) < c(i + i/ c y+c\b\tfdtE 2 (t) (415) 

when c = a/b. It is interesting to note that the above conditions have basically the same 
form as (|4.11|) , apart from the factor e ( 1+1 / c )i+ c (~ 0.62 in the MSSM.) 
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5 Preliminary applications and comments 



The aim of this section is to illustrate briefly through some examples the possible use of 
the integrated forms and make contact with the existing studies and approximations. It 
is, however, obviously not meant to be exhaustive nor refined from the phenomenological 
point of view (for instance no threshold effects or higher loop effects are included), as this 
would deserve a separate analysis by itself. Let us also keep in mind that all the functions 
Ei which enter the solutions are analytically known in terms of the initial gauge coupling 
values as can been seen from eq. ( |1.1| ) and its solutions. 



5.1 Landau pole free bounds: 

The necessary and sufficient bounds found in the previous section are exact generalization 
of the one in [@] initially derived in the regime of large yti^> yb)- The necessary bounds 
( [4.11 ) which would be also sufficient in the limit a — > ( where they become identical to 
( f4.14| |4.15| ) ) give a first estimate of the allowed values for y t , y b at low energy, restricting 
them to a rectangle 



V? SSM {EW) < y t , yt ISSM {EW) < y b (5.16) 

where y t and y b depend on the gauge couplings and are easily determined numerically, 
for instance in terms of a grand unified scale, the value of the gauge couplings at that 
scale, and some low energy electroweak scale. 

For instance, one has roughly [neglecting the g\ and g 2 couplings], 



r -_ 7 9j{EW) 

Vt > b ~ 18i_ f 4fflw/9 (5 - 17) 



We stress that the necessary bounds in eq.( 5.16| ) involve no approximation whatsoever. 



One can readily turn them into constraints on the tan (3 parameter at some electroweak 
scale: 




The bars indicate that the masses and gauge couplings are running quantities at the 
chosen low energy scale. 

Going now to the improved bounds ( |4.12|) one can reduce further the allowed range for 



the Yukawa couplings. These bounds do not allow in their general form an easy analytic 
determination of the allowed regions. To get a feeling about these regions, let us illustrate 
the case in two different approximations: 

i) Ei = E 2 = E, i.e. neglect the difference f 1 — f 2 = gf/24ir 2 
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ii) assume yb(to), Vt(to) sufficiently small for a first order expansion to be legitimate 



i) In this case one integral can be performed exactly in fl4.12|) , leading to 



Vt{U) < 
Vb(t ) < 



lfr(to)(l-f) 



l-(l-\b\y b (t )$dhE(t 2 ;t )y-*/ b 

yt(t )(i-f) 

l-(l-\b\y t (t )$dt 2 E(t 2 ;t )y- a / b 



(5.19) 



The domain defined by ( |5.19| ) lies within the rectangle ( |5.16| ) and is controlled by the 
relative strength of a and b. For instance, the delimiting curves start off at the points 
(y b , 0) and (0,y t ) with first derivatives equal to — ^ ( — ^ in the MSSM and — | in the 
SM) , to be compared with in the rectangular approximation. Moreover, the exact fixed 
line solution y t = yb leads to the constraint 



yt{ta) = Vbito) < 



E(t ) 



\a + b\fpdtE(t) 



(5.20) 



to be compared with the rectangular approximation (|5.16|) where \a + b\ is replaced 
by (a 14% effect in the MSSM and a 25% effect in the SM.) These considerations give 
a qualitative guideline of the reduction of the allowed domain. 

ii) In this case one gets a linearized approximation of the domain [keeping though the 
effect of the U{\)y coupling], in the following form, 



' \b\(f E i y yt (t ) + \a\(JE 1 JE 2 )y b (t ) < J E x 

(5.21) 

k \b\(J E 2 fy b (t ) + \a\(JE 2 J E x )y t {t Q ) < J E, 



where the E^s here are normalized to Ei(to), J ... = J* T ° ...dt\ and 
f...J... = #...dt } $...dt 2 . 

In this approximation the necessary domain is delimitated by two straight lines with 
scale dependent slopes. Again, one can translate these conditions into bounds on tan/? 
at some effective electroweak scale. 



5.2 yt — yb ~ 93 approximation, fixed points and quasi- fixed line: 

In the approximation where Ei(t) = E 2 (t) = E(t), and assuming that the initial values 
yt,y® are small enough so that one iteration in the form of eqs. (|2.12|j2.13|) is a good 
approximation for the Eij(t), an easy integration yields 
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yt(t) 
y b (t) 



-by° b JE) a / b 

y° b E(t) 



1 



1 + 

1 



(5.22) 
(5.23) 



provided of course all other Yukawas are put to zero. Note that at this level of 
approximations the above solutions depend just on one integral, namely Jq E. If one goes 
further and neglects the SUl(2) gauge coupling then Jq E can be computed explicitly, and 
one obtains, 



Pt(X) = p° t 



X 



1 



Pb{X) = Pbj 



1 + apl(X - l)j [1 + (1 + apl{X - 1) )i-c - 1) 

X 1 



1 + ap° t (X - l)j |l + (1 + ap° t (X - 1) )i- - 1) 



(5.24) 
(5.25) 



where now we use the reduced variables 

Vt{t) 



Pt 



pb 



y b (t) 



and 



93 



(5.26) 



(5.27) 



(Note that a = 18/7, c = a/b = 1/6 in the MSSM and a = 9/2, c = 1/3 in the SM.) 
The approximate solutions ( |5.24| , |5.25| ) allow one to retrieve the well-known infrared fixed 
points in the (p t , pb) plane, fl7|| . For instance, in the MSSM, the fixed point (pt = Pb = 
0) ( respectively (p t = 0, pb = y^) ) is obtained by looking at the limiting behaviour of 
pt{X) when p° — > (resp. of pb(X) when p° — > 0. On the other hand, the IR (attractive) 
fixed point (p t = \, Pb = |) is obtained by expanding Pt(X) and Pb(X) simultaneously for 
small p° and p°. It should come as no surprise that those exact fixed points are obtained 
only in this limit since the solutions ( |5.24j , |5.25| ) become exact only in this limit, [see also 
a related comment at the end of this section.] 

The exact fixed line p t = Pb is trivially obtained from ( |5.24| , |5.25|) , i.e. starting 
from p° = p° the two reduced Yukawas remain equal at any other scale. Perhaps more 
interesting is to ask whether one can determine analytically the other exact or effective 
IR fixed lines. An IR attractive effective fixed line of the form 



Pt + Pb = - 
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(5.26 



was found in |T§ for the MSSM. 

Starting from eqs.( [5.24] , |5.25|) one can actually improve on this effective fixed line in 
the following way. For small Pt,p b one obtains the integral line 



a{l - c)(p° b - p° t )p tPb + p°(l - a(p° b + cp° t ))p b - p°(l - a(p° t + cp° b ))p t = 



(5.29) 



If we require this integral line to go through the exact fixed point pt = Pb = a h +c ) 
then a fixed line is obtained in the form 

Pt + Pb 2(l + c)« (5.30) 



i 



Ptpb 

to first order in an expansion around the fixed point pt = Pb — a (x+ c ) ■ 
That is 



= 6 (5.31) 

Ptpb 

in the case of the MSSM, which constitutes an improved effective IR fixed line beyond 
Eq. (|5.28|) . One can even reasonably expect Eq. (|5.31|) , and more generally Eq. (|5.30|) , to 



be exact in the regime under consideration. Indeed, for instance in the case of the SM, 
the effective fixed line ( |5.30|) reads 



= 12 (5.32) 

Ptpb 

On the other hand, an exact IR fixed line is known in this case (see |fPT|] ) and is the sum 



of two terms, one of which coincides precisely with ( 5.32 ), the other being vanishing to 



first order in the deviation, 5, around the fixed point, that is for p t = 1/6 + 5, pb = 1/6 — 5. 

Let us end this section by noting further possible applications of the integrated forms. 
As we mentioned before, if one starts from ( |5.24| , |5.25| ), one retrieves the exact fixed point 
Pt = Pb = l/( a (l + c)) only in the region of small Pt,Pb- Moreover, in the deep infrared 
region (X — > oo) one obtains from Eqs.( |5.24 |5.25| ) p t = Pb = (1 — c)/a irrespective of 



the initial values p®,p b - The IR attraction to this point is of course an artifact of the 
approximate solutions. In fact, one can resum exactly the integrated forms ( |2.12| , |2.13j ) 



in the limit X — ► oo and obtain pt = Pb = 1/( Q! (1 + c )) as the attractive IR fixed point, 
independently of the initial values Pt,p b - 

Other regimes can be also looked at (for instance p° p b or p b p°) for which 
approximate analytical expressions for the integrated forms can be obtained up to three 
iterations, thus improving on ( |5.24| , |5.25| ). This allows to tackle the form of the fixed line 



in such regimes. We do not dwell further on these aspects here. 
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5.3 constraints in the y t — — y T and all gauge couplings case: 

In this section we illustrate the use of the solutions in the t — b — r system only to derive 
inequalities which correlate the three fermion running masses in terms of initial values 
for the Yukawa couplings. Starting from equations ( [2.16| - |2.18|) and using the fact that 



[1 + \ot\ifi fui)'^ > 1 for all expressions of that form appearing in Q2.16| , [2.17D one writes 
immediately the following (necessary) inequalities: 



m t 



< Y t °(2V2G F )- 1/2 ^E' 1 sm(3 (5.33) 



< Y?{2V2G F )- 1/2 JEiCosf3 (5.34) 



where i = 2(6), 3(r) and recalling that the y[s are squares of the Yukawa couplings Y{. 
Y® denote the values of these couplings at some high energy scale and the bar denotes 
running quantities at the electroweak scale. The E[s are as defined in eqs.([2.9|, |2.18|) . 



These inequalities give boundary conditions on the initial values of the Yukawa couplings, 
necessary to retrieve the correct "physical" fermion masses. For instance one immediately 
sees from (|5.34j ) that large tan/3 necessitates large initial values for the bottom and r 
Yukawa couplings. 

Moreover, relying systematically on the fact that (1 + (aji/j / Wj)'^' > 1 for i = 1,2,3, 
one can derive the following optimal rigorous inequalities for more involved combinations 



(5.35) 
(5.36) 
(5.37) 



yl 8 


> 


(y° b ) 18 




yhf 


(ymy° T 


) 35 EfE'i 


yhr 


> 


(y?ny°r 


fEfEi 


yl 2 


(y° b ) 12 


Ef 


yt 


> 


{yl) 4 


El 


yfyr 


(y?) 21 y° T 


EfE, 



which can be readily translated into inequalities involving the running quark masses, 
tan/3 and the three gauge couplings (all taken at the electroweak scale), as well as the 
values of the three Yukawa couplings at some initial scale 



fnE 18 ^ (Y, ) 18 El (2V2G 



- /^OW^flVW / _„ „:„3 a 17 a ^°- 0i5 



,10 



mT 3 ^ 35 ~ (Y 4 °) 3 (Y r ) 35 / e 3 E? 5 sin 3 {3 cos 17 f3 



WW (y t 0) 2 (y T 0)3 El tan2 ^ 

~ (y0)i2 V^3^|^(2V2G F ) (5.39) 

^ 4 > (Y b r El (2V2G F f 

~ m°) 21 n° /^T^ tan 3 /? sin 18 /? K ' ] 



m t 21 m 7 
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These inequalities express general necessary conditions which delineate the physically 
allowed regions for the initial values of the three Yukawa couplings, i.e. consistent with the 
values of the physical top, bottom and r masses, prior to any model assumption^. [Note 
also that (|5.33| , |5.34j) are already contained in ( |5.38| - |5.40|) .] Finally, had we neglected 
the r Yukawa coupling, the necessary inequalities involving the top and bottom running 
masses would have read 



m t > y t ° y^taji/3 r- 5/2 
~ (Y b ) e El cos 5 /3 l V f) 

m? ~ (Y t °) e Ef tan/? sin 5 /T F> 



(5.41) 
(5.42) 



The different origin and meaning of these inequalities as compared to ( |5.18| ) should be 
clear. 



5.4 A numerical illustration 

Even though the general form of the exact solutions is not directly exploitable analytically, 
truncated iterations provide very good approximations which can be furthermore very well 
controlled using the convergence criteria we derived. For instance truncating at the first 
iteration, i.e. approximating E\2 and E 2 % in Eqs.( (p.l2 , [2.13|) by the explicit forms 



E 12 {t) 
En(t) 



Ex{t) 



[l - byl !% E 2 {t')dt'Y' h 
E 2 {t) 

:i-by?foEi(t>)dtr /b 



(5.43) 
(5.44) 



and plugging them back into Eqs. (|2.10| , |2.11|) , one gets a simple analytical solution 



In table 1, a comparison is made for this solution with the Runge-Kutta method, showing 
an excellent agreement of less than 1% accuracy for any small, moderate or large values 
of tan/3. 

Similar approximations can be obtained for the top-bottom-r case, at least if the initial 
Yukawa couplings verify the sufficient convergence criterion of section 3.2. We will not 
dwell on further possible applications in the present paper. 



6 Conclusion 

We have written down integrated forms for the running of the Yukawa couplings in the 
case of two and three Yukawa fermions, which are easily generalizable to any number of 

4 In a more refined treatment, one should of course correct for the difference between the running and 
the pole masses 
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tan (3 


y0 
1 b 


Y? 


Y b (t), truncated 


Y b (t), R.-K. 


Y t (t), truncated 


Y t (t), R.-K. 


2 


0.0387453 


1.13007 


0.0145059 


0.0145050 


0.775788 


0.775974 


10 


0.174138 


1.01581 


0.0630978 


0.0631052 


0.54263 


0.542743 


50 


0.866544 


1.01097 


0.435682 


0.439526 


0.585453 


0.590258 



Table 1: Numerical comparison between the exact one-loop solution (truncated to the 
first iteration) and the Runge-Kutta RG evolution. The evolution is over 10 orders of 
magnitude starting from the initial Yukawa coupling values shown in the table 



such fermions. These forms are exact solutions for the one-loop renormalization group 
equations, valid for virtually any gauge theory with a Yukawa sector. The most important 
feature of such forms is that they allow for a rigorous determination of convergence criteria 
as well as exact conditions for avoiding Landau-like poles of the Yukawa couplings. In 
the case of the MSSM, such criteria lead to approximate analytical solutions in the top- 
bottom system, with very good numerical accuracy ( < 1%) for any value of tan f3. Similar 
criteria were obtained for the top-bottom-r system, which lead to controllable analytical 
approximations. In this context we gave some preliminary applications for Landau pole 
bounds in the top-bottom system, commented on some infra-red fixed points and lines, 
and gave optimal necessary constraints on the Yukawa couplings both in the top-bottom 
and top-bottom-r systems. 

In view of the increasing phenomenological interest for the large tan (3 scenario, such 
solutions should prove useful in determining the exact structure of the running of the 
remaining parameters of the MSSM using for instance the method developed in |Tj|, and 
possible implications on the structure of the (stable) infra-red fixed points [pOR . Very 
recently, the authors of reference |2TJ have addressed similar issues, starting though from 
approximate solutions. 



Acknowledgment: We are indebted to J.-L. Kneur for providing us with the nu- 
merical illustration presented in section 5.4 and thank him as well as C. Le Mouel for 
discussions. This work has been performed partly in the context of GDR-Supersymetrie 
where preliminary results were published in . 
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A: Necessary and sufficient conditions for non singu- 
lar evolutions 



A:.l A necessary condition 

A necessary condition not to meet a singularity in running up from an initial energy scale 
t = to to a given high energy scale t = T < t Q is easy to establish. Indeed, if yb,t(t) are 
free from singularities in the interval [T, to], then E^it; t ) are necessarily positive for any 
t in this interval, as can be seen from eqs.( [4.6|) and the fact that b < and y b ,t{t) > 0. It 
then follows from eqs.([4.7|) that 



£12 (Mo) > Etfato) and E 21 (t,t ) > E 2 {t,t ) (A.l) 



since the denominators in eqs. ( f4.7|) are always smaller than one (recall that a/b > and 
t < t ). From the above considerations one gets immediately the inequalities 



\b\y t {t ) f ° dt'E l2 {t'; t ) < 1 - |% t (* ) f° dt'E 1 (t'; t ) (A.2) 
Jt Jt 



1 - \b\y b (t ) f dt'E 21 (t';t ) < 1 - \b\y b (t ) f° dt'E 2 {t'-t 
Jt Jt 



(A.3) 



for any t in the interval [T, to]. Again, from eq.( [4.6| ), the lefthand side of eq. ( |A.2|) should 
remain positive for any t in the interval [T, to], for yt{t) being free from singularities there. 
In particular t = T gives the most significant condition, whence 



1 - \b\y t (t ) £° dt'E^t'-to) > 



(A.4) 



which is the necessary condition given in eq.( |4.11| ). One obtains similarly the second 
inequality in eq. (|4.11| ). 

On the other hand, from ( |A.3|) and the positivity of its left-hand side one gets 



Ei 2 (t; t c 



Exit; t c 



E\(t; t 



(1 - by b (t ) Jl E 21 (t>; t Q )dtr /b ~ (1 - by b (t ) g E 2 (t>; to)dt'Y' h (A _ 5) 



valid since E% and a/b are both positive. One can thus repeat the same proof which 
lead to (|A.4|) with Ei(t;t ) replaced by 



E\(t; t c 



(l-by b (to)Jt E 2 (t>;t )dt>)*/t> 

in eq. (|A.l| ) to get the first inequality in eq. ( [4.121 ), and the second in a similar way. The 
infinite tower of inequalities (|4.11[ - [4. 13|) is obtained recursively in the same way. 
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A:. 2 A sufficient condition 

Similarly to what was done in section 3, eqs.( |4.7|) define a mapping A in the form 



e'u(*) 

4(0 



1 

[I - \b\y° b It' e 2 (t>)e 21 (t>)dt>)°/» 
1 

;i-|%H t0 ei(f)e 12 (f)df) a /6 



where 




(A.6) 



c (*) = |4^4 (A.7) 

e;(t) = ^(t;to) (A.8) 

2& = M*o) (A.9) 

Again, we collect the ey(t)'s in a vector 



E(t) = ei2 (A.10) 



and consider the range T < t < to where to corresponds to some low energy scale at 
which initial values for y t , y b are chosen, and T a high energy scale (typically a GUT scale) 
up to which we require the Yukawa couplings to have a non singular behaviour. We also 
define a norm similar to eq.(|3.6|) 



|| E ||= max{ sup | ei2(t) |, sup | e2i(t) |} (A. 11) 

T<t<t Q T<t<t 

To determine the conditions we are looking for to avoid singularities in the range 
[T, t ] it will actually suffice to ask when does the mapping defined in eq. (|A.6[) become a 



contraction. Parts of the the proof will resemble that of section 3.1. However, in contrast 
to the latter case where the mapping defined in eq. (|3.4j ) could not have singularities as 



long as Ui(t),U2(t) > 0, in the present case one has to make sure that the mapping A 
keeps ejj(t) within a finite interval 1 < e^t) < R. 

For a given R, let us thus denote by X R the set of all vectors E(t) such that 1 < 
eij{t) < R for any t in the interval [T, to]. We look for conditions on the values of R, y®, yl 
such that 

• i) The mapping A sends any element of Xr in Xr, (so that the ey(t)'s remain in 
the interval [1, R] after an arbitrary number of iterations of A.) 
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• ii) A is a strictly contracting mapping in Xr, that is 

|| A{E X ) - A{E 2 ) \\< K R || E x - E 2 || (A.12) 

with some Kr < 1. 
Condition i) means that 

4j{t) = il'\b\yl t ^e^)e 3l {t')dtri b ~ R (A ' 13) 
On the other hand, one finds from eij(t') < R that 

e ' < ( a 14^ 

13 ~ {l-\b\yl t R^e 3 (t')dtr/ b 1 ' ] 

provided that 1 - \b\y% t R //° ej(t')dt' is a positive numberL]. 
In view of eqs. ( |A.13| , |A.14j ) a sufficient condition to obtain i) is 



< R 



ll-\b\y?RS}°e 1 (t>)dt')°/« 

and 

1 



(l-\b\y° b Rf* e 2 {t')dt'Y' h 
which translates easily into 



1 1 

yi < 



< R 



R R^ 1 ^ 
1 1 



where c = a/b and 



(A.15) 



yi = \b\ij!j\i(t)dt (A.16) 



At this level R is still an arbitrary number. However, the optimal situation would be 
to chose it such that the upper bound ( |A.15| ) be the largest possible. This would be the 



5 This condition will, however, turn out to be already contained in the sufficient condition we are 
looking for. 
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case for R = (1 + l/ c ) c > but one still has to check whether this value is compatible with 
the second requirement ii) which we turn to now. 
Using the inequality (|3.8|), one gets from eq. ( |A.6| ) 



e'n{t)- e' l2 {t) | < 



c\b\y$fT dt'e 2 (t') \ e 21 {?) - e 21 {?) 



[1 - \b\y^inax{S^dfe 2 (t')e 21 {t')J^dt'e 2 {f)e 21 {t')}Y- 



< c\b\y° b ^dt'e 2 (t') ^ | 

" [l-|6|y 6 it:^ e 2 (f)^'] 1 + c 11 



(A.17) 



valid for any i in the interval [T, to] and E, E belonging to Xr. A similar inequality holds 
obviously for | e 21 (t) — e' 2l {t) |, and one finally gets 



A(E) — A(E) \\< K R || E-E 



(A.18) 



with 



Kr = max{ 



cy2 



cyi 



{1-Ry 2 y+^ (1-R yi ] 



l+c 



} 



(A.19) 



where y\,y 2 are defined in eq. ( |A.16| ). It is now easy to check that when condition 
( |A.15|) is satisfied with strict inequalities, one gets 



K R <1 



(A.20) 



even for the value of R quoted before, Rq — (1 + l/c) c , which maximizes the bounds 
in eq.( [A.15[ ). Since eqs.( |A.l8[ |A.20[ ) mean that the mapping is indeed contracting, one 
concludes that the sufficient conditions for i) given in eq. ( |A.15| ) with maximal bounds, 
i.e. 



2/2 < 



1 



c(l + l/c) 1+c 

1 

c(l + l/c) 1+c 



(A.21) 



imply also ii). It follows that when (|A.21|) (equivalently ( [4.14 , 4.15 ) ) are satisfied, a 
unique, regular, solution for the eqs.fl4.7p exists in X Ro . The regularity of y t , yj, as given 
by eq. ( [4.6|) is then an immediate consequence. 
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fB.r 



B: Approximate solutions for y t ^> ^ 

To first order in y® one finds for y t 



v (t) = y ° El(t) h + ay ° b I* E ^' ] dt> 

where 



= /' Exti'W (B.2) 

JO 

and where the solution for is given in eq. fl2.7|) . 

C: Exact integrated forms for an arbitrary number 
of Yukawa couplings 

Under the restriction of flavour- conserving Yukawa couplings, and assuming that the Higgs 
fields sit in representations such that the renormalization group equations for the Yukawa 
couplings can be cast in the following form at the one-loop level |Lf| : 

j t Vi = Vi(fi(t) + OijVj) (CI) 

where i, j count the fermion fields, and yi denotes the square of the i th Yukawa coupling. 
Then the exact solution for each reads: 

Vi(t) = : V>1 ft (C2) 
1 - a-uVi Jo u i 

where the u^s are given by the implicit equations 

u .(t) = E ^ (C3) 



and^(t) = e/o^(^'. 
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